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An interesting dynamical model for dark energy which does not require extremely light scalar fields 
such as quintessence, and at the same time explains the (near-) coincidence between the neutrino 
and dark energy densities is the model of dark energy coupled to mass varying neutrinos (MaVaNs). 
Despite the attractions of this model, we show that, generically, this model contains a catastrophic 
instability which occurs when neutrinos become non-relativistic. As a result of this instability, as 
neutrinos become non-relativistic, they condense into neutrino nuggets which redshift away similar 
to cold dark matter, and thus cease to act as dark energy. Any stable MaVaNs dark energy model 
is extremely contrived, and is virtually indistinguishable from a cosmological constant. 



I. INTRODUCTION 

During the past decade, it has become clear that dif- 
ferent cosmological observations, such as the dimming 
of distant supernovae la [11121, anisotropics in the cos- 
mic microwave background J3] , and the large scale struc- 
ture of the universe (e.g., jj|) cannot be explained with 
a cosmological model that contains only ordinary (bary- 
onic+dark) matter evolving according to Einstein's the- 
ory of general relativity. The most popular solution is 
to introduce an extra component with negative pressure, 
the so-called dark energy (e.g., Q)- 

While the simplest candidate for dark energy is a cos- 
mological constant (or vacuum energy) the extremely 
small size of the cosmological constant is difficult to ex- 
plain unless one is willing to invoke a selection effect (e.g., 
|6j). Furthermore there is an additional strange coinci- 
dence, that the density of photons, neutrinos, baryons 
and dark matter and dark energy density all seem to be- 
come comparable at roughly the same time in the history 
of the universe even though they scale very differently 
with redshift. 

Dynamical models of dark energy (e.g., quintessence 
0, k-essence H, ||) may be built to try to explain the 
coincidence. These models require an extremely light 
scalar field (with mass < 10~ 33 eV). An interesting al- 
ternative, which simultaneously alleviates the need for 
an extremely light scalar field, and explains the coinci- 
dence between the neutrino and dark energy densities 
has been recently proposed in the context of a theory of 
mass varying neutrinos (MaVaNs) [ToL ITT| (also see [l2|). 
In this model, dark energy is modulated by the density 
of neutrinos, which can simultaneously cause the evolu- 
tion of both dark energy density and neutrino masses 
as a function of cosmic time. The time evolution how- 
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ever is not caused by a slowly rolling light scalar, but 
by a scalar (the so-called acceleron, A), which is heavy 
(relative to the Hubble expansion rate) but has a min- 
imum that varies slowly as a function of the density of 
neutrinos. The model may also have interesting impli- 
cations for neutrino oscillation experiments |10l [l3j and 
solar neutrino observations [lj, [l5j . 

In this work, we investigate the evolution of pertur- 
bations in the MaVaNs model, which has not yet been 
consistently studied in the literature. The key differ- 
ence between perturbations in MaVaNs and most other 
dark energy candidates is that (unlike quintessence or k- 
essence) all the dynamical properties of (non-relativistic) 
MaVaNs are set by the local neutrino density. In par- 
ticular, pressure is a local function of neutrino density, 
implying that the hydrodynamic perturbations are adi- 
abatic. Subsequently, the speed of sound for adiabatic 
perturbations at scales smaller than the Hubble radius 
is set by the instantaneous time derivatives of pressure 
and density . We will show that this implies that hy- 
drodynamic perturbation are unstable on all macroscopic 
scales once the neutrinos become non-relativistic. 

The paper is organized as follows: in Sec^ we sum- 
marize properties of the homogeneous MaVaNs model 
for dark energy. Sec lIIII presents hydrodynamic and ki- 
netic theory arguments for why perturbations of non- 
relativistic MaVaNs have an imaginary sound speed, and 
are thus unstable. In particular we will show that, when 
neutrinos are non-relativistic, free streaming is unable 
to stabilize perturbations on all scales with wavelengths 
larger than m^ 1 where m_4 is the mass of the acceleron 
field. We will then move on to speculate about the 
end result of the instability. In Sec II VI we argue that 
the outcome of instability is inevitably a multi-phase 
medium with most of neutrinos in a dense phase, which 
we call neutrino nuggets, and dilute as pressureless mat- 
ter. SecE] presents a thermodynamic description of the 
phase transition and draws the analogy with the liquid- 
gas phase equilibrium. Finally, Sec lVIl considers poten- 
tial stable MaVaNs (or adiabatic dark energy) models, 
and Sec IVIII concludes the paper. One thing will become 
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clear: After neutrinos become non-relativistic, the model 
has a serious microscopic instability. As a result of this 
instability, the fluid will cease to act as dark energy and 
will not be able to drive the acceleration of the universe. 



II. PRELIMINARIES OF THE DARK ENERGY 
WITH MASS- VARYING NEUTRINOS 

Within the MaVaNs theory of dark energy, the local 
energy density associated with dark energy is the sum of 
the energy densities of neutrinos and a light a scalar field, 
A, the so-called acceleron field, which also modulates the 
neutrino mass: 



V = n v m v (A) + V (A), 



(1) 



where the neutrinos are assumed to be non-relativistic. 

In the limit n v 3> ni A , where tua is the mass of the 
acceleron field, the acceleron field responds to the average 
density of neutrinos, and relaxes at the minimum of the 
potential V (the actual relaxation of A at the minimum, 
of course, depends on its evolution through the cosmic 
history). In this limit, both m v and Vq become functions 
of the neutrino density, n„ , through 

dV ( dV Q \dm v dV 

if dm v /dA ^ 0. 

Assuming that the kinetic energy in the acceleron 
field is negligible [2(j , the equation of state for the neu- 
trino/ acceleron fluid is: 



Pressure 



w 



V (A) 



-1 



n v m v (A) 



(3) 



Density V V 

where we used the fact that non-relativistic neutrinos 
have negligible pressure. We note that w ~ — 1 (i.e. a 
cosmological constant), as long as the energy density in 
neutrinos is small enough compared to the energy of the 
acceleron field. 

Fardon et al. go on to propose a specific model for 
the mass varying neutrinos: 

o 

m v (A) 



V (A) = A 4 ln[l 



M{AY 
■M(A)/n]+C, 



(4) 
(5) 



where it is assumed that A4(A)/(i ^> 1, and A ~ 10 _3 eV 
characterizes the energy scale of dark energy and C is just 
a constant needed to set the minimum of the potential to 
zero (i.e. no cosmological constant). Now, using Eq.|J2J, 
we find: 

V ^ K Hn(^)+C^ niJ = — , (6) 



where mo 
state: 



mf r /p, and Eq.© gives the equation of 



w = — 1 



(7) 



Based on current observational constraints |1 + w\ < 
0.2 (e.g., Q), and thus m > 10 2 to„. Therefore, we 
restrict the rest of our analysis (with the exception of 
Sec. [Vj to the tuq ^> m v regime. 

For the model to be sensitive to the mean density of 
neutrinos, there needs to be many neutrinos within the 
Compton wavelength of A, i.e. N v ~ n u /m A ^> 1. The 
mass of A is given by, 



fd 2 V\ 



A4_ 

2 



7/7 



di 



OA 



We thus have 



A 2 



dA 



-3 



(8) 



(9) 



Observational bounds restrict the coupling (dm u /dA) < 
10 -6 01 which puts a lower bound on N u around 10 24 , 
justifying the assumption that A field is only sensitive to 
to the mean neutrino number density. 

For non- vanishing dV/dm v , Eq.© can also be satisfied 
if dm v jdA = for a value of A or m v . This can happen 
if Vo(*4) reaches a minimum, which subsequently puts 
a maximum on m v < m max . m max corresponds to the 
mass of neutrinos in vacuum or when ?!„ < A 4 /m max . 
In particular, the macroscopic dynamics of the original 
MaVaNs model jll] is quantified in terms of A, toq and 



m 



V (m v ) = A In H + C, m v < m 



(10) 



Again, C = —A 4 ln(l + mo/m max ) is added such that the 
minimum of the potential is at zero. 



III. EMERGENCE OF AN IMAGINARY SPEED 
OF SOUND 

A. Hydrodynamic picture 

In the conventional quintessence models 0, similar 
to inflationary models, the scalar field is slowly rolling 
at the present epoch, and therefore its effective mass is 
smaller than the Hubble expansion rate. In contrast, 
the acceleron field sits at the instantaneous minimum of 
its potential, and the cosmic expansion only modulates 
this minimum through changes in the local neutrino den- 
sity. The mass of the acceleron field can be much larger 
than the Hubble expansion rate. However, by the same 
token, the coherence length of acceleron, m^ 1 , is much 
smaller than the present Hubble length, and thus, un- 
like quintessence, the perturbations on sub- Hubble scales 
> vn~^ are adiabatic, i.e. obey the same equation of state 
as the homogeneous universe. Therefore, the speed of 
sound, c 2 , for these perturbations is simply given by |16| : 



P 
9 



2 -i Wp + Wf) 
C-Z = — = : = W 



d\w, 



3H{l + w) d\nn v 



(11) 
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where we used the continuity equation p = — 3H (l + w)p, 
and in the last step, we used Eq.Q. 

We can already see that for a vanishing or small enough 
w, c 2 has the same sign as w, which is negative. Clearly, 
this is a catastrophic instability in this theory, as the sub- 
Hubble perturbations on scales larger than m^ 1 (H < 
k < m_4) would grow as exp(fc|c s |f) if c 2 < 0. 

In particular, for the model introduced in ^lj, Eq.® 
gives: 



P 



A 4 



const. 



P 
9 



(12) 



B. Kinetic theory picture 

One may wonder if, similar to the gravitational insta- 
bility of collisionless cold dark matter (or massive neu- 
trinos with constant mass), our instability is stabilized 
at small scales due to free streaming of neutrinos. In 
this section we study the nature of this instability in 
the kinetic theory picture, and show that the imaginary 
speed of sound is indeed scale-independent in any non- 
rclativistic MaVaNs models, making the instability most 
severe at the microscopic scales (~ m^ 1 ). 

The action of a test particle with a dynamical mass 
m„(x, 77), in a flat FRW universe, is given by: 



S = — J m v (x, r/)ds = 
where 



m l/ (x, rj)a(r))\/l — u 2 dr], 
(13) 



u • u 



dx <ix 

drj drj ' 



(14) 



and a and r\ are the scale factor and conformal time re- 
spectively. Thus, the single particle Lagrangian, Momen- 
tum, and Energy take the standard form 



C v = —m v a\f\ 
dC v m„au 



du 



2 

u ■ p - C v 

where 7 = (1 — v 



= m„7au, 



2\-l/2 



(15) 

(16) 

(17) 
(18) 



while we have dropped the x, r\ dependence for brevity. 

The evolution of p is given by the Euler-Lagrange 
equation: 



drj 



<9x 



(19) 



where V denotes the gradient with respect to the comov- 
ing coordinates x, and we have neglected gravity. 

Now the Boltzmann equation for the evolution of the 
phase space density of neutrinos, /(x, p), takes the form: 



df 
drj 



u • V/ — 07 l Vm v 



d£ 
dp 



0. 



(20) 



Let us also generalize Eq. (J5J) for relativistic neutrinos. 
Using the fact that neutrino momentum remains constant 
in a homogeneous background, the relaxation of the A 
field at the minimum of its potential leads to 



&V 
OA 



=/(x,p)d 3 p + 



y/p 2 + ml ' drriv 
which for dm v /dA 7^ (m v < m max ) yields 



~dA 



= 0, 



drrir. 



A 4 



(21) 



(22) 



for the Logarithmic potential. 

In the absence of perturbations, neutrinos become non- 
relativistic as they are free-streaming. Therefore, as a 
result of momentum conservation, they are frozen into 
the relativistic Fermi-Dirac phase space density: 



/(P) 



2(2tt 



,-3 



2(2tt)- 



exp[|p|/r ,„] + l exp \p ph /T„] + 1 ' 



(23) 



where p p h = |p|/a is the physical momentum, while 
T v = T a u /a is the neutrino kinematic temperature, 
which matches the thermodynamic temperature in the 
relativistic regime. The factor of 2 in Eq. (|23|l accounts 
for both neutrinos and anti-neutrinos. 

Let us study the linear perturbations of the Boltz- 
mann equation (Eq. I20fl . in the sub-Hubble regime, 
i.e. for time/length scales much shorter than the Hubble 
time/length. We will show that the instability occurs on 
microscopic scales so this is perfectly adequate. In this 
regime, without loss of generality, we can assume a = 1. 
For plane- wave linear perturbations: 



Sf = A(p) exp [i(k ■ x - uji])] , 
8m v — £ exp [i(k ■ x — wr/)] , 



Eq.lO gives: 



df 



cA(p) = (k.u)A( P )- 7 - 1 ^k.^J S, 
while 8VJ dA = constraint (Eq. I21f> yields 



r 1 A(p)d 3 P - 



AS = 0, 



(24) 
(25) 



(26) 



(27) 



where 



d 2 V 
dm 2 , 



d 2 v 

dm 2 



in. 



(p 2 - 
(p 2 



il) 3 ' 2 



f(p)d 3 p 



01JL 



(28) 



We first notice that that there is no preferred scale 
in Ea. (|2*S|l . as only the combination uj/k (= c s ; speed of 
sound) , appears in the equation. This is contrary to what 
happens for gravity where there is a scale, the Jean's 
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scale, below which random motions stabilize perturba- 
tions. The reason for the difference is that for scales 
k < nij[ the acceleron field adjusts to the local density 
rather than satisfy a Poisson like equation. 

Taking the first term on the right hand side to the left, 
dividing by u> — k • u, and substituting from Eq. 1(27(1 
yields: 



A( P ) 



k-tT)/7-V)A(p')d 3 P 



ku)A 



(29) 



We can eliminate the integral of the unknown amplitude 
A(p) through multiplying both sides by 7 _1 (p), and inte- 
grating over d 3 p, which give sthe characteristic equation 
for c, 



A = 



d 3 p 



l 2 (p)(c s - k ■ u) 



dp 



(30) 



In the limit of u 2 <C \c^\, 1, we can expand the argu- 
ment of the integral in Ea. (|30|) . in powers of u, and ignore 
0(u A ) terms. The angular parts of the integrals can be 
taken using the spherical symmetry of /(p), which yields: 



(Ai:p 2 dp) 

3 3 
ic! ~ 2 



P 



im v c 2 

2 



_P_ 

rn„ 



\ 4 

01 JL 



luC 2 



1 + (c s 



5/2)4 
mf, 



O 



_P_ 

m„ 



(31) 



where, in the last step, we used integration by parts. 
Combining Eqs. l(2"5|) and (jBTj) yields: 



l + (c 2 s + c s 



mf, 



O 



P 

m„ 



(32) 

Finally, following Eq. 122(1 , the relation between n v and 
77i„ is also corrected at finite temperature: 



1 



(P 2 



2m 2 



dVp 



(33) 



which can be substituted in Eq. (|32|) to give: 

2 ^dVo fd 2 v \ 

c„ = m.. - — ^ 



V) 



(34) 



Note that in the limit of p <ti m u , this result is equiva- 
lent to Eq. 1)11(1 for non-relativistic neutrinos, as riu ~ 
—dVo/dm^ in this limit, so that 



_i dV 



dm v \ dm 2 



d In) 



d In n v 



= -1, 



(35) 



where the last equality is for the logarithmic potential. 

Moreover, we should point out that the sign of the 
sound speed is set by the sign of the second derivative of 



V(A) at the extremum, d 2 V/dA 2 



d 2 V Q d 2 V d 2 V (dm v 
dm 2 



dml dA 2 V dA 



We have 



for p 2 < ml, (36) 



so that 



m v m\ 



dm v 
~0A 



(37) 



Therefore, any realistic MaVaNs scenario with m\ > 
becomes unstable to hydrodynamic perturbations in the 
non-relativistic regime. 

In particular, for the Fermi-Dirac neutrino phase space 
density (Eq. E3J), we have {p 2 )/T 2 = 15C(5)/C(3) ~ 12.9, 
and assuming the Logarithmic model of Sec. ITU Vn ~ 
— lam^ + const., we end up with 



= -1 + 51.8 — 



m. 



O ( — 

m,. 



(38) 



Thus, we see that for the logarithmic potential, the model 
is unstable if T v < to„/7.2. Using Eq. ©, it is easy to 
see that the unstable regime corresponds to T v < 1.2A. 



IV. THE OUTCOME OF THE INSTABILITY: 
NEUTRINO NUGGETS 

As we argued above, there is no preferred macroscopic 
(> ttCa ) scale in the linear problem. In particular, the 
free-streaming of neutrinos acts uniformly on all macro- 
scopic scales, and may only stabilize the perturbations of 
relativistic neutrinos. Therefore, the inevitable outcome 
of the instability is the formation of non- linear structures 
in neutrino density or neutrino nuggets. As the instabil- 
ity is the fastest at smallest scales we expect nuggets to 
form instantly at scales larger than m~^ , after the onset 
of instability. The non-linear nuggets may then slowly 
merge to form larger nuggets. However, one may also en- 
vision a scenario in which non-relativistic neutrinos are 
continuously heated back up to relativistic temperatures, 
as a result of instability, and at the expense of reducing 
the vacuum energy of the A field. We investigate this 
possibility below: 



A. Is it possible to have a homogeneous and 
marginally stable phase? 

Assuming that this scenario works, neutrinos would be 
constantly heated up at microscopic scales through scat- 
tering off non-linear neutrino nuggets. As the instability, 
and hence scattering occurs at microscopic scales, neutri- 
nos will have a microscopic mean free path, which implies 
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that evolution should be adiabatic on macroscopic scales. 
The assumption of marginal stability then implies: 

p 

c 2 = — = ^ P = const, (39) 
P 

which is necessary to maintain homogeneity and marginal 
instability at the same time. However, using energy con- 
servation, p + 3H(p + P) = 0, we can see that p and P 
can only have a ACDM-like behavior: 

P = A+-,, (40) 
a 6 

P=-A, (41) 

where A and B are integration constants. In particular, 
we require that 

u 2 B 

p + P = n„m,„( 1 (l + —)) = —, (42) 

O (X 

where we used the fact that P — n^m^^u 2 ) /3 for the 
neutrino gas. However, using the fact that n„ oc a -3 , we 
find 

, . u 2 , B B , . 

w„(7(l H )) = = const. => m u < . (43) 

3 n , M no <v 

The fact that m„ reaches a maximum, implies that either 
Vo(m y ) (or Vo(-4)) reaches a minimum, or (7 -1 ) oc n" 1 ex 
a 3 (eq. 1221 • The latter, however, can only last for less 
than a Hubble time, as 7 ~ 1 at the onset of instability, 
while (7 _1 ) < 1. 

Therefore, we see that the only way to maintain ho- 
mogeneity of marginally stable neutrino gas is for the 
acceleron to settle at its true minimum, following which 
the neutrinos are decoupled and become non-relativistic, 
while the evolution follows the standard ACDM scenario. 



B. Production of Neutrino Nuggets: Qualitative 
description of a multi-phase outcome 

In Sec. IIV Al we argued that the instability cannot be 
followed by a homogeneous (single phase) state of the 
neutrino-acceleron fluid, unless the acceleron settles at 
its true minimum (i.e. vacuum). If this does not happen, 
the outcome can only be a multi-phase medium. 

As result of the instability, the neutrino fluid should 
fragment into condensed neutrino nuggets (or the liquid 
/interactive phase) , within which, neutrinos maintain rel- 
ativistic energies to have a stable density distribution. 
The remaining neutrinos remain in a tenuous gas phase, 
and should remain relativistic to maintain their stabil- 
ity. Moreover during the instability there is a possibility 
of converting part of the energy into A particles. The 
A particles are not trapped inside the nuggets and their 
contribution to the energy density diminishes with the 
Hubble expansion. 



As neutrinos in the "gas phase" continuously lose en- 
ergy to Hubble expansion, and nuggets, due to their small 
volume fraction, cannot significantly affect the gas phase, 
most of the neutrinos should accrete into the nuggets. In 
fact, as we show in Sec. [V] within the assumption of 
thermal equilibrium, this process instantly exhausts the 
neutrinos outside the nuggets, unless they can decouple 
(i.e., A settles at its true vacuum) within the gas phase. 

The fraction of neutrinos in the gas phase is set by the 
equilibrium of evaporation rate from and accretion rate 
onto the surface of the nuggets. However, the large den- 
sity contrast between inside and outside of the nuggets 
implies that only neutrinos with a large Lorentz factor 
(7 > m out /m m ~ n in /n out ) can escape the nuggets. As 
(7) inside nuggets is limited by a few times its initial 
value at the onset of the instability (where (7) was ~ 1), 
the equilibrium can only be maintained when a small 
fraction of neutrinos remain outside the nuggets. 

Thus we conclude that the most probable outcome of 
the instability is the formation of dense non-linear struc- 
tures of neutrinos, surrounded by practically empty space 
where the acceleron has settled to its minimum. All this 
is happening at microscopic scales. As far as the expan- 
sion of the universe is concerned, the fluid energy density 
will redshift as matter and thus cannot drive the accel- 
eration of the Universe. 



C. Production of A particles 

Another possible outcome of the instability is produc- 
tion of A particles. In fact, one may speculate that, since 
the characteristic time and length scale for the instability 
of non-relativistic MaVaNs is m^ 1 , a significant fraction 
of neutrino/ acceleron energy may end up in A particles 
(or accelerons). However, we should note that, at the 
onset of the instability, the neutrinos are still marginally 
relativistic, implying that the characteristic time for the 
instability may be significantly longer. 

Since the speed of sound (and thus the instability) is 
modulated by the Hubble expansion, it is reasonable to 
approximate c 2 as: 

c 2 = -a(HAt) + O(HAt) 2 , (44) 

close to the onset of the instability, where a ~ 1, is con- 
stant of order unity. The maximum growth occurs at the 
scale of Ti^ 1 , and is given by 

A oc exp(|c s |m^Ai) = exp(<xff 1/2 At 3/2 m A ). (45) 

Therefore, the characteristic growth time for the insta- 
bility is given by 

At ms ~ (^) 1/3 m ^ > m A \ (46) 

since, by construction, acceleron is much heavier than 
Hubble scale in MaVaNs theories (see Sec. [HJ|. This 
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implies that the production of A particles cannot be a 
significant energy sink, at least during the initial stage of 
instability (or nugget formation). 

Subsequent mergers of relativistic m~^ 1 -sized nuggets 
may potentially lead to A particle production. However, 
the above argument can also be used to show that the 
velocities of produced nuggets, v nug , can never be rela- 
tivistic: 



< 



At, 



< 1. 



(47) 



NUGGET-GAS PHASE TRANSITION: 
THERMODYNAMIC PICTURE 



THE 



In Sec. IIV Bl we presented a qualitative argument for 
why most of neutrinos must condense into nuggets as 
they become non-relativistic. In this section, we present 
a more rigorous picture of the associated phase transition 
within the assumption of thermodynamic equilibrium. 

However, we should first note that, in the homogeneous 
phase, the thermodynamic approximation is clearly in- 
valid, as the cross-section for neutrino interaction is 
small, and thus, by construction, thermal relaxation 
takes much longer than the Hubble time However, 
in the relativistic regime, the distribution remains ther- 
mal as the expansion does not change the original shape 
of the thermal Fcrmi-Dirac distribution. After neutrinos 
become non-relativistic, the instability develops within 
microscopic time and length scales, and thus, it seems 
plausible to assume that the system reaches some sort of 
thermal equilibrium within microscopic times. 

Assuming thermodynamic equilibrium, the phase 
space density of neutrinos is given by the Fermi-Dirac 
distribution: 



dn,, 



d 3 p 



exp 



T 



i — l 



1 



d 3 p 



(48) 

where the additional factor of 2 accounts for both neu- 
trinos and anti-neutrinos, p, is the chemical potential, 
which can be non-zero if the neutrino number is con- 
served, while T stands for the thermodynamic tempera- 
ture. 

Similar to previous sections, we assume that the ac- 
celeron field, A, is relaxed at its minimum. Therefore, 
the system is subject to the constraint condition of Eq. 
(|21|l . which in conjunction with the condition: 



d 3 p 
(2tt) 3 



/(p), 



(49) 



set the equation of state of the system, as they fix all ther- 
modynamic variables (e.g., [i and m u ) for given neutrino 
number density, n u , and temperature, T. In particular, 
the pressure of the system is given by: 



P = ^n^m^^u 2 ) - V Q (m v ). 



(50) 



Sh 

w 

0) — (C 
u 



i 1 1 1 1 1 r 




_L 



-2 2 4 6 

-log n(A 3 ) 

FIG. 1: Pressure- Volume phase diagram for the Logaithmic 
MaVaNs model with mo = 10 3 A and m ma x — > oo. The solid 
curves are isothermal contours with T = 2 A, for integer ts, 
where the highest temperature is T — 1024A. The isothermal 
contours become degenerate below T ~ 0.5A. The thermal 
average of neutrino Lorentz factor, (7}, is equal to 2 on the 
dotted curve, separating the relativistic and non-relativistic 
regimes. The dashed curve is the adiabat with /1 — > as 
T — > 00, which approximates the cosmic history of neutrinos 
in the absence of inhomogeneities. 



ID 

a. 



-2 - 




2 4 

-log n(A 3 ) 

FIG. 2: The same as Fig. gj for m max = 10 2 A. 



7 



Fig. shows the pressure as a function of specific vol- 
ume (= n^ 1 ) for the system described by these equations 
for the logarithmic potential of Sec. |n]with mo = 10 3 A 
(Eq. I1(J[1 . The solid curves show different tempera- 
tures of T„ = 2 e A for integer Ps, where the highest 
plotted temperature is 1024 A. The curves become de- 
generate below T v ~ 0.5 A, as the degeneracy pressure 
dominates the thermal pressure. The dotted curve cor- 
responds to (7) = 2, which separates the relativistic and 
non-relativistic regimes. 

A homogeneous medium becomes unstable if its pres- 
sure decreases with increasing density. In Sec. IIII Bl 
we saw that this happened when neutrinos become non- 
relativistic, which is also clearly seen in Fig. (JIJ, where 
dP/dn changes sign (i.e. isothermal curves reach a mini- 
mum) at the boundary of relativistic and non-relativistic 
regimes. 

The dashed curve in Fig. JIJ is the adiabat with /i = 
as T — > 00, which approximates the cosmic history of Ma- 
VaNs for the logarithmic potential. We again see that 
the pressure reaches a minimum (rendering the system 
unstable), when the neutrinos become non-relativistic, 
which also happens close to where the adiabat crosses 
T = A isothermal curve, as demonstrated in Sec. IIII Bl 
Therefore, as the density of neutrinos decreases mono- 
tonically due to cosmic expansion, they become unsta- 
ble and fragment into nuggets at n v ~ A 3 , as no sta- 
ble phase is available at lower densities. This can also 
be seen in the asymptotic (large n v ) limit, as the ki- 
netic pressure is simply given by n^T, which is sub- 
dominant to the acceleron pressure for m v ^> m , as 

— Vo = — A 4 ln(l + mo/m„) oc — nj 2 , using the mass- 
density relation in the relevant range. Therefore, dP/dn v 
remains negative at smaller densities, and so the system 
does not have any other stable phase, justifying the qual- 
itative picture of Sec. IIVI 

A more interesting case of a finite maximum neutrino 
mass, with m max = 10 2 A is shown in Fig. J2J). The 
curves are the same as in Fig. (Q. Here, we see that for 
densities lower than a critical density (which depends on 
temperature), the neutrino mass is fixed at its maximum, 
implying that the acceleron field relaxes at its true vac- 
uum. Therefore, we end up with a gas of free-streaming 
neutrinos with the acceleron settled at the minimum. As 
a result, the system shows two stable phases of interactive 
(or liquid) mass-varying neutrinos and non-interactive 
(or gas) neutrinos with m u — m max . The two phases 
can co-exist only if they have the same pressure, temper- 
ature and chemical potential. 

Fig. Q is a blow-up of the T = 16A isothermal con- 
tour from Fig. J3J). Points A and C on the isothermal 
contour are liquid and gas phases which can coexist at 
this temperature. A pictorial way of finding the coexis- 
tence states is via Gibbs construction. Introducing e and 
s, as the energy and entropy per unit particle, the ther- 
modynamic and extensivity relations subsequently take 




0.5 1 1.5 2 2.5 3 
-log n(A 3 ) 

FIG. 3: A blow-up of Fig. iffi. which shows the Gibbs con- 
struction for T — 16A. Points A and C are interactive and 
non-interactive (liquid and gas) phases, which can coexist as 
they have the same pressure, temperature and chemical po- 
tential. This requires areas Sab and Sbc to be equal (after 
including the logarithmic metric in the plot), as explained in 
the text. The dotted line corresponds to (7) = 2, while the 
lower solid line shows the T = degenerate curve. 

the form 

de = Tds- Pdn' 1 , (51) 
e = Ts- Pnr 1 + fi, (52) 

which can be combined to give |17| 

d f i = rT x dP - sdT. (53) 

Integrating over the isothermal contour ABC, the condi- 
tion of equal chemical potentials for the coexisting phases 
A and C takes the form: 

- = / n~ 1 dP = Sab - Sbc = 0. (54) 
J ABC 

Therefore, the area enclosed in regions Sab and Sbc 
must be the same (even though it may not look this way 
due to the logarithmic scale in Fig. ©. 

A. Onset of phase transition: bubble nucleation vs. 
linear instability 

Even though the neutrino-acceleron fluid does not be- 
come unstable to linear perturbations until it becomes 
non-relativistic, it can become unstable to perturbations 
of finite amplitude, otherwise known as bubble nucleation, 
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much earlier. We will argue however that this is not rel- 
evant in the cosmological context being studied. 

Consider a bubble of true vacuum with size larger than 
m^ 1 in the relativistic interactive (liquid) phase. If the 
pressure of the ambient medium is smaller than the vac- 
uum pressure inside the bubble (which is set to zero in 
our treatment), the bubble rapidly expands until it is 
in pressure equilibrium with the ambient medium. This 
can in principle be achieved through contraction of the 
medium as a result of bubble nucleation. 

However, the onset of bubble nucleation depends on 
the statistics of seed bubbles present in the system. In 
particular, the smaller the value of m^, the larger will be 
the minimum size of the bubble that can trigger a phase 
transition. The bounds on m_4, which come from the 
requirement of no thermal acceleron production in early 
universe , imply more that 10 24 neutrinos in a volume 
of mT? , making it extremely unlikely for seed (vacuum) 
bubbles to exist in a Hubble volume. Cosmological fluc- 
tuations of neutrino fluid are also damped as a result of 
free-streaming, and thus cannot seed bubble nucleation. 
Therefore, the phase transition is unlikely to start un- 
til the pressure reaches its minimum, where small linear 
perturbations can seed the instability. 

The outcome of the instability, which develops on mi- 
croscopic time scales, can be uniquely obtained from con- 
servation of total energy and particle densities. For ex- 
ample, for the case of m max — ► oo where pressure and 
density vanish in the vacuum, the final density of frag- 
mented nuggets is roughly at the intersection of /i = 
adiabat (dashed curve in Fig. ^) with the P = axis, 
which is approximately where T v ~ Afh^mo/A)] 1 / 4 . 



B. The fate of nugget-gas multi-phase state 

Finally, let us consider the fate of the system of coex- 
isting neutrino nuggets and gas. Assuming that mo > 
Tn max > A[21|, the area of the region Sab (see Fig. |3J) 
is dominated by the non-relativistic regime and is given 
by: 

S. i u - ! n-HP ~ m max + O ( A, . (55) 



AB 



mo 



Neglecting the degeneracy pressure (which can be eas- 
ily justified for m rnax > A), the area of region Sbc is 
given by: 



Sbc = 



n Tdn 



-T\n(n gas m max / A 4 ). 



(56) 



Therefore, the Gibbs construction Sab — Sbc yields 



L gas ' ""max 



A exp 



T / 



for T < A < 



« A 4 



(57) 



As we anticipated in Sec. IIV Bl as long as the phase 
transition occurs far enough from the true minimum of 



Vo(m v ), the density of neutrinos in the gas phase (where 
A relaxes at its true vacuum) is exponentially suppressed. 
The subsequent evolution follows from requiring energy 
conservation: 



d(n„a 3 T) = —n u a 3 m max d 



<Tda J 



1 



ln(n ffas m r . 



C M 4 ) 
r da 3 



\u{n gas m max / A 4 ) ' 
which can be approximately solved to give: 
T 1 



(58) 



'gas 



In (7 



gas 1 

1 



,/A 4 ) 



m„ 



c /A + 31n(a/a A ) : 



(59) 



where n v denotes neutrino cosmic mean density, and a\ 
is the cosmological scale factor at phase transition (where 
T ~ A). Therefore, both the fraction of neutrinos in the 
gas phase, as well as the equilibrium temperature, remain 
essentially unchanged their cosmic history. 



VI. DISCUSSION 
A. How about a light neutrino mass eigenstate? 

We have shown that unless the acceleron field sets the 
mass of a relativistic particle, hydrodynamic perturba- 
tions are unstable. The Big Bang model predicts a kine- 
matic temperature of T v ~ 0.7 Tq mr ^ 1-7 x 10 -4 eV for 
relic cosmological neutrinos (e.g., [13). Current exper- 
imental limits on neutrino oscillations require the most 
massive neutrino eigenstate to have m„ > 10 -2 eV fl9| . 
Therefore, at least one of the neutrino eigenstates for 
cosmological neutrinos is non-relativistic. 

One may imagine a scenario where the interaction with 
A is suppressed for the non-relativistic (including the 
most massive) eignestates, and thus the acceleron/dark 
energy dynamics is modulated by the density of a rel- 
ativistic neutrino eigenstate. Even this seems difficult. 
For the logarithmic potential, the stability of this sce- 
nario requires T v > 1.2A (Sec. IIIIB|I . which is already 
constrained, up to a logarithimc factor, by the current 
dark energy density (Q DE — 0.7): 

(2.3 x 10~ 3 eV) 4 = p DE ~ A 4 In (m /m u ) 
< 0.5 T 4 In (m /m„) ~ (1.4 x 10~ 4 eV) 4 In (m /m„)(j30) 



requiring 



\n(m /m u ) > 7 x 10 4 



(61) 



which does not seem realistic. 

Even if we do not restrict Vb(m„) to a logarithmic form, 
assuming a Fermi-Dirac distribution (Eq. 12311 . we have 



0.6 T 4 

Pde 



10- 



(62) 
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for any relativistic MaVaNs theory, implying a very slow 
evolution that is indistinguishable from the ACDM sce- 
nario. One may hope to be able to test this model by 
measuring neutrino oscillations (or structure formation) 
at high rcdshifts, as neutrino mass scales as the inverse of 
the temperature of neutrino background in such models 
(m„ oc T^ 1 oc (1 + z)~ x ). 

B. Adiabatic dark energy and cosmological 
constant 

In Sec. IIV Al we showed that an adiabatic dark en- 
ergy model with c s = behaves exactly like a ACDM 
cosmology. The argument can be easily generalized for 
any stable dark energy model (c 2 . > 0) which satisfies the 
null energy condition (p DE + P DE > 0) in the following 
way. Combining the stability condition and conservation 
of energy we have 

P P 

J2 de * de . n (ao\ 

c s = -■ = -TF77 T~p \ > °> ( 63 ) 

P de 3H(p DE +P DE ) 

which can be integrated to give 

/•CO 

P DE (t Q ) - 3 / H(t)c 2 s (t) [p DE (t) + P DE {€)} dt 

J to 

<P DE (t o )<0. (64) 

Since dark energy has negative pressure at the present 
time (P DE (to) < 0), Eq. (fo4"j) guarantees that the pres- 
sure remains negative and finite, even when all parti- 
cle densities approach zero at future infinity. There- 
fore a stable adiabatic dark energy model will necessarily 
asymptote to a cosmological constant/vacuum energy at 
future infinity (or it goes unstable). 

VII. CONCLUSIONS 

We have demonstrated that the model of dark energy 
with mass- varying neutrinos f MaVaNs) [Tol Hl| is subject 



to a linear instability with an imaginary speed of sound. 
The onset of instability is around the time that neutrinos 
become non-relativistic, and is due to the fact that, unlike 
quintessence perturbations, the fluctuations in the non- 
relativistic neutrino density evolve adiabatically. The ac- 
celeron field mediates a strong attractive force (stronger 
than gravity) between neutrinos. There are more than 
10 24 neutrinos inside the Compton wavelength of the A 
field. When neutrinos become non-relativistic, their ran- 
dom motions are not sufficient to stop collapse and the 
medium becomes unstable. 

The question of the outcome of the instability, an 
inherently non-linear process, is a more difficult prob- 
lem to address. We present qualitative arguments for 
why the outcome of the instability should be a multi- 
phase medium, where most of neutrinos end up in dense 
nuggets (interactive/liquid phase) with nearly constant 
density, while a small fraction may remain outside in a 
decoupled gas phase. We then confirm these arguments 
quantitatively within the assumption of thermodynamic 
equilibrium, which is plausible following the phase tran- 
sition. Both nuggets, and free neutrinos (which are expo- 
nentially suppressed) dilute similar to cold dark matter 
particles, with negligible pressure, leaving a vacuum en- 
ergy /cosmological constant as the only possible source of 
the observed cosmic acceleration after the onset of insta- 
bility in the MaVaNs models. 
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